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^! Abstract 

We perform numerical simulations of the Hubbard model using the projector 

Quantum Monte Carlo method. A novel approach for finite size scaling is 

Q ' discussed. We obtain evidence in favor of d-wave superconductivity in the 

repulsive Hubbard model. For U = 4, T c is roughly estimated as T c ~ 30 K. 
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After the discovery of high-temperature superconductivity (HTSC) the two-dimensional 
Hubbard model JT|, M has been proposed as a model for a theoretical explanation of the phe- 
nomena. Indeed it has been shown, that the Hubbard model exhibits the normal conducting 
and magnetic properties of HTSC ||. 

It is now widely accepted, that HTSC show d-wave symmetry of the superconducting or- 
der parameter 0,0. According to our simulations || , and recent work || the repulsive 
Hubbard model also favors d-wave symmetry. But the question of (d-wave) superconduc- 
tivity in the Hubbard model has been discussed controversially ||, @, |§. 

We proposed the tt '-Hubbard model as the suited model for numerical simulations ||, 
0. It exhibits a Van Hove singularity away from half filling 0. Furthermore with the 
tuning of the next nearest neighbor hopping t' we are in the position to circumvent the 
majority of numerical difficulties in the simulation. The tt'-Hubbard model is described by 
the Hamiltonian 

H H Ub = - t j:i: cU> - E E c U, CT + f/ E ^n id . (i) 

<i,j> CT <Cij> CT i 

Here c\ a creates an electron with spin a on site i, n ita is the corresponding number operator 
and U is the on-site Coulomb interaction. The sum <i,j> (^C i, j ^>) runs over the pairs 
of (next) nearest neighbors. 

Our simulations are performed with the projector quantum Monte Carlo method 
(PQMC) §, 13, U. In which the ground state 

|*„> = ^- 6H \^t) (2) 

of the Hamiltonian 7i is projected from a testfunction Y&t) with a normalization constant 
M and with the projection parameter 9. Details of the method are described in [ jT^] . 

To provide evidence for superconductivity we follow the standard concept of off diagonal 
long range order (ODLRO) |13| . Therefore we study the vertex correlation function 



C V d (r) = - E E 00* (cl T ci +5i c i+r+n c J+rT ) + yEE 969*<T*(i, r)Cf(i + 5, r + 5') (3) 

i 5,5' i 5,5' 

with the phase factors gr^, grj ± 1 to model the d-wave symmetry. As shown in || and in 
further detail in M the d-wave correlations are positive for larger distances \r\ and level off 
to a "plateau". Other superconducting symmetries (in particular s-wave) fluctuate around 
zero. This results has recently been supported by f|. Our current simulations reach the 
same conclusion for the pure Hubbard model (£' = 0). 

The question of superconductivity can be only answered by finite size scaling. In the case 
of weak or intermediate interaction |14]] the behavior of the correlation function is dominated 
by the shell structure of the system. Considering the average vertex correlation function 

C^jE^Ar) (4) 

with the number of lattice points L = L 2 X standard 1/L X scaling for instance seems to 
provide clear evidence against superconductivity |16|], ||. In this paper we argue that 
this conclusion is too simplified. 

In this context we introduce a BCS-reduced Hubbard model, the J-model, with the 
same kinetic Hamiltonian as the tt '-Hubbard model and an interaction favoring cooper 
pairs with d-wave symmetry. In momentum space the Hamiltonian of this model is given 
by the Hamiltonian 

h bcs = E E £ k4ck + 7 E /fc/?4,T ct -M c -M c P,T (5) 

k cr=±l k,p 

with the single particle energies = — 2t( cos(k x ) + cos(k y )) — 4t f (cos(k x ) ■ cos(k y )) and 
the form factor = cos(k x ) — cos(k y ) for modeling the d-wave interaction and = 1 
for the s-wave interaction. Superconductivity has been rigorously proven for this type of 
BCS-reduced Hubbard models [DJ . 



Figure [l] shows the 1 / L x scaling of for the J-model. For weaker interaction we would 
reach the same conclusion, the absence of superconductivity, as other authors in the case of 
the Hubbard model. Only the case J = —0.25 would be superconducting. Considering the 



susceptibility (Xd = LC%) again for weaker interaction the divergence of \d ls ambiguous. 
This results of figure [I] have been obtained with the stochastic diagonalization (SD) ||18|| , 



TJ, [0. Details are published in 0, 0. 



To avoid these problems with the standard finite size scaling ansatz we propose a novel 
approach J7J. Finite size scaling is provided by comparation of the Hubbard model to the J- 
model. Before we carry out this comparation we have to circumvent a further complication 
in the Hubbard model. The values of the correlations functions C^(r) are extremely large 
for smaller distances |r| ||, and therefore susceptible to the fluctuations in the numerical 
simulations. Indeed these fluctuations for smaller distances exceed the "plateau" value of 
the C^(r). As only the long range behavior is of interest for superconductivity we restrict 
the average vertex correlation function 

0^ = 1 £ C» (6) 

Li c r 

\r | > | r c 

to the distances \r\ > \r c \. In equation || |r c | is a critical distance and L c is the number of 
points with |r| > |r c | for r e 1, . . . L. Typically we choose |r c | = 1.9. 

We now make the assumption, that the same finite size scaling behavior of the supercon- 
ducting correlation functions C^' p is equal for both models. Accordingly we have for each 
interaction U and system size L an effective J e ff of the J-model. 

We determine J e ff in the following way: For the system parameters L, (n) (filling) and 
t' and the interaction U we calculate C v d ' p for the Hubbard model with PQMC For the same 
set of parameters we tune J using the SD method to obtain the same value C^ p in the 
J-model. This J is our J e ff- 

A first test was carried out for the negative (attractive) Hubbard model, which is com- 
monly believed to be superconducting (s-wave symmetry). Results in table | show a unique 
J e ff for system sizes 6 x 6 to 12 x 12 and small deviations at 4 x 4. |r c | = 1.9 was chosen 
for this and all following cases. It should be mentioned, that the choice of \r c \ is not critical 
for the qualitative behavior of J e ff- 

In table [lj we return to the repulsive Hubbard model. For U = 2 and t' = —0.22 we 



again obtain as in the attractive case a unique J e ff for system sizes 6 x 6 to 12 x 12. We 
notice a decrease for 16 x 16. The same effect occurs for the case U — 1. But for t' — (the 
pure Hubbard model, U = 2) we find agreement up to 16 x 16 (table |T1[). 

This behavior is explained by the existence of different finite size gaps. For the small 
gaps in the t' = —0.22 case the simulations are only valid for relatively large projection 
parameters 9, which exceed our numerical possibilities. In figure § we show the increase of 
C v d ' p with various 9. In contrast the upper curve shows the leveling off in the t' = case for 
a still moderate 9. This is caused by the relatively large finite size gap. 

Therefore we conclude, that the deviation for 16 x 16 in the if = —0.22 system is due to 
insufficiently large 9 in the simulation. The system does not reach the ground state properly. 
Larger 9 are outside of the reach of methods. At this point we would like to mention that 
energy measurements are still rather insensitive to 9 compared to the vertex correlations 



12| . This is a rather important point as agreement in energy measurements was often used 
as evidence for the validity of a certain numerical method. 

Considering again table |TJ we conclude that our simulations show clear evidence for the 
existence of d-wave superconductivity in the Hubbard model. 

The simulations had to be restricted to values U < 2 and system sizes up to 16 x 16 
because of the convergence problems of the PQMC method outside this parameter regime. 
This is clearly indicated by a dramatic break down of the average sign (table H). Figure |3| 
shows the regime of "safe" simulations below the shaded areas. 

The effective interaction J e ff leads to a superconducting T c in the BCS-model |ZT]. The 



BCS-T C has to be considered as at least a rough estimate and it does not include fluctuations 



in the two-dimensional system. Simulations by Schneider et al. |22| suggest that for the 
range of our interactions the deviation of the BCS-T C and the real T c is rather small. Our 
simulations clearly suggest a considerable difference between T* (forming of pairs) and T c 
as described by Schneider et al. J/J, fl23| . 

For the J eff « 0.07 in table [TJ [U = 2, f = -0.22) we only obtain a very low T c « 1 K. 
But in a 6 x 6 system we are able to calculate U — 4 with a sufficient large 9. For t' = —0.22 



we find J e // ~ 0.2. This effective interaction leads T c pa 30 K. For t' = we obtain a T c pa 5 K. 
J e jj « 0.2 agrees very well with a recently published 12 x 12 lattice ||. The difference of T c 
is explained by the closeness of the Van Hove singularity in the t' = —0.22 case 0. Larger 
values of C v d ' p for U — 6 and [/ = 8 as suggested by 4 x 4 exact diagonalization results may 
lead to a dramatic increase of T c . But we do not want base this decision on 4 x 4 lattice 
sizes. The conclusion of || that the Hubbard model does not exhibit superconductivity for 
larger U and larger L is not valid. The errorbars of Cd{r) are about ten times larger than 
the value of C v d predicted by our J-model simulations. 

In conclusion we provide clear evidence for the existence of d-wave superconductivity 
in the Hubbard model. For U = 4 we obtain a T c pa 30 K. Therefore the single band 
Hubbard model has to be considered as a serious candidate for the explanation of high T c 
superconductivity. 

We want to thank T. Husslein, D.M. Newns, H. De Raedt, T. Schneider, J.M. Singer 
and E. Stoll for the helpful discussions and ideas. This work was supported by the Deutsche 
Forschungs Gemeinschaft (DFG). The Leibnitz Rechenzentrum (Munich) grants us a gener- 
ous amount of CPU time on the IBM SP2 parallel computer. 
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FIGURES 
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FIG. 1. Finite size scaling of the averaged (left) and the cumulated (right) vertex correlation 
function. The ground state of the BCS-reduced Hubbard model with d-wave interaction was 
calculated with the stochastic diagonalization. The filling of table 2 and t' = was used. 
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FIG. 2. The ^-scaling is plotted for the 16 x 16 system with the filling = nj_ = 105, the 
interaction U = 1 and if = —0.22 (runs: t' = —0.22) and for the 16 x 16 system with the filling 
n-f = nj_ = 109, the interaction U = 2 and if = (runs: t' = 0). The average sign (sign) is in both 
cases about one and 6/m = r = 0.125. (a) 9 verses ground state energy (b) verses the averaged 



vertex correlations function with \r r \ = 1.9 C 
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FIG. 3. Limits of the PQMC for the tt'-Hubbard model with the filling (n) » 0.8 and f = / 
t' = -0.22 
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TABLES 



L 


72t = 72 l 


u 


t' 


c v 


J *> f f 


4x4 


5 


-0.5 


0.0 


0.00196(4) 


-0.38 


6x6 


13 


-0.5 


0.0 


0.00089(1) 


-0.31 


8x8 


25 


-0.5 


0.0 


0.00056(2) 


-0.30 


10 x 10 


41 


-0.5 


0.0 


0.00038(1) 


-0.30 


12 x 12 


61 


-0.5 


0.0 


0.00028(1) 


-0.30 


4x4 


5 


-1 


0.0 


0.0045(1) 


-0.76 


6x6 


13 


-1 


0.0 


0.00218(3) 


-0.64 


8x8 


25 


-1 


0.0 


0.00146(2) 


-0.63 


10 x 10 


41 


-1 


0.0 


0.00105(1) 


-0.64 


12 x 12 


61 


-1 


0.0 


0.00079(1) 


-0.64 



TABLE I. Effective interaction J e of the BCS-reduced Hubbard model with s-wave interaction. 
The PQMC has used 6 = 8 and m = 64. The statistical errors of the last digit are given in the 
brackets. 
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L 


"T = n i 


U 


f 


(sign) 


r< v 'P 


J eff 


6x6 


13 


2 


0.0 


1.000 


0.00134(2) 


-0.077 


8x8 


25 


2 


0.0 


0.999 


0.00099(2) 


-0.075 


10 x 10 


41 


2 


0.0 


0.995 


0.00073(2) 


-0.073 


12 x 12 


61 


2 


0.0 


0.987 


0.00056(2) 


-0.071 


16 x 16 


109 


2 


0.0 


1.000 


0.000278(6) 


-0.078 


6x6 


13 


1 


-0.22 


1.000 


0.00053(2) 


-0.025 


8x8 


25 


1 


-0.22 


1.000 


0.00037(1) 


-0.025 


10 x 10 


41 


1 


-0.22 


1.000 


0.000259(4) 


-0.023 


12 x 12 


61 


1 


-0.22 


1.000 


0.000193(4) 


-0.022 


16 x 16 


105 


1 


-0.22 


1.000 


0.000104(4) 


-0.016 


6x6 


13 


2 


-0.22 


0.999 


0.00134(2) 


-0.061 


8x8 


25 


2 


-0.22 


0.995 


0.00107(2) 


-0.069 


10 x 10 


41 


2 


-0.22 


0.973 


0.00080(2) 


-0.068 


12 x 12 


61 


2 


-0.22 


0.876 


0.00063(1) 


-0.067 


16 x 16 


105 


2 


-0.22 


0.567 


0.00023(4) 


-0.026 



TABLE II. Effective interaction J e of the BCS-reduced Hubbard model with d-wave interac- 
tion. The PQMC has used = 8 and m = 64. The statistical errors of the last digit are given in 
the brackets. 
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